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In graphene, the pseudospin and the valley flavor arise as new types of 
quantum degrees of freedom due to the honeycomb lattice comprising two 
sublattices (A and B) and two inequivalent Dirac points (K and K') in the 
Brillouin zone, respectively. Unique electronic properties of graphene result 
in striking phenomena such as Klein tunnelling, [U, |2] Veselago lens,|3j and 
valley-polarized currents. [4 — 7J Here, we investigate the roles of the pseudospin 
and the valley in electron scattering at graphene edges and show that they 
are strongly correlated with charge density modulations of short-wavelength 
oscillations and slowly-decaying beat patterns. Theoretical analyses using 
nearest-neighbor tight-binding methods and first-principles density-functional 
theory (DFT) calculations agree well with our experimental data from the 
scanning tunneling microscopy (STM). We believe that this study will lead to 
useful application of graphene to "valleytronics" and "pseudospintronics". 

A variety of edge properties [5HTT] of graphene and graphene nanoribbons[12j have been 
investigated and interference images using the STM were also reported before. [13] However, 
electron scattering behaviours at graphene edges have not been well understood yet. A con- 
ventional metal with a terrace and a step can be thought as accomodating a two-dimensional 
(2D) free electron gas with a hard wall edge and the standing wave formed at the edge can 
be analytically solved. This behaviour was directly observed at the steps of Au(lll) and 
Cu(lll) surfaces by scanning tunneling microscopy (STM). [T4lfI(3] Now, a question arises as 
to whether the graphene edge has a similar standing wave pattern to conventional metals, 
and how two sublattices and two inequivalent valleys in graphene affects the scattering and 
the standing wave formation. We show below that the interference pattern at the graphene 
edge is quantitatively understood in terms of intra- and intervalley scattering processes at 
the graphene edge and that local defects at the edge considerably change the interference 
pattern in certain cases. 

Due to the crystal momentum conservation along the edge, available backscattering chan- 
nels are limited to the Bloch states of the same wavevector component in the edge direction 
(k y ) as the incident wave. The two distinct Fermi circles around K and K' valleys in 
the doped graphene are folded back onto the rectangular Brillouin zone (equivalent to the 
hexagonal first Brillouin zone) in the armchair and the zigzag edge cases, as depicted in 
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Figs, la and lb, respectively. For a given wavevector k y , there are two intersecting points 
for each Fermi circle, corresponding to left and right propagating states. Figures lc and 
Id show two possible (intravalley and intervalley) backscattering channels at the armchair 
edge whereas only one intravalley scattering channel is allowed at the zigzag edge. [17] In 
the nearest-neighbor tight-binding method, the scattering wave at graphene edges can be 
uniquely determined by boundary matching conditions. As shown in the left panels of Fig. 
2, we choose three representative model structures [T8] for the armchair edge termination 
that have been observed in the transmission electron microscopy image [TU] as well as an- 
other structure of the zigzag edge termination. The resulting reflection probabilities and 
electronic density profile are calculated with appropriate hopping parameters from the first- 
principles Wannier function analysis (See Supplementary Information) and the Fermi energy 
(Ep) is chosen to be 0.3 eV above the Dirac point (Ed)- Note that the epitaxial graphene 
on the SiC(OOOl) surface is n-doped and typically Ep — Ed ~ 0.3 - 0.4 eV.[T9] 

Considering the crystal momentum conservation at the armchair edge, intervalley and 
intravalley scatterings can in principle be allowed as mentioned above (See Fig. lc). How- 
ever, when the boundary matching conditions are applied (See Supplementary Information), 
the actual scattering at the armchair edge is an entirely intervalley process for all incident 
angles. [20] and the pseudospin of the incident wave is identical to that of the scattering wave 
as can be inferred from Fig. la. In other words, the pseudospin is invariant throughout 
the scattering process at the ideal armchair edge. (In the pseudospin formalism, A and B 
sub lattice components of the wavefunctions are described as up and down "spins".) By inte- 
grating all scattering waves on the Fermi surface, we construct laterally averaged electronic 
density profile in the real space (in the right panel of Fig. 2). For the armchair edge, the 
electronic density profile can be written as 

f\ i#» + ^(x)i 2 + 1^>) + ^r\*)\ 2 ^ 

= Jjv | exp(i(— K — k x )x) — exjp(i(K + k x )x))\ 2 d9 
+ /A I exp(i(— K + k x )x) — exp(i(K — k x )x))\ 2 d6 

oc 1 — cos(2K x) J (2k R x) , (1) 

where x is the distance from the edge, ^(jcn*^ O^) is the incident (scattering) planewave on 
the K(K')- valley Fermi surface with the angle 9, K the length of TK in fc-space, k R the 

3 



radius of the Fermi circle and Jo the zero-order Bessel function. 

Very interestingly, the armchair edge has intervalley scattering which is not shared by the 
2D free electron gas. In the electronic density profile near the armchair edge, a short-period 
oscillation and a slowly decaying beat pattern are observed. First, a short-period oscillation 
comes from the intervalley scattering, and its wavelength is approximately 1.85 A, i.e., 3a/4 
where a is the length of graphene lattice vector (0.246 nm). Therefore, this phenomenon 
is a unique property of materials with multiple Fermi circles. Second, the slowly-decaying 
oscillatory behaviour has the same origin as the one observed on the terrace of conventional 
metal. pH] Its characteristic period is a few nanometres. The Bessel-function-type envelope 
pattern is in general caused by quantum interference between incident and scattered electrons 
in 2D metals at the straight edge termination. Finally, a beat is an interference between 
two waves of slightly different wave vectors. As depicted in Fig. 1(c), | — A^i | (K to K' 
scattering) is slightly different from \k 3 — [K ! to K scattering), and these two intervalley 
scattering events result in the beat. Since the beat is originated from the finite radius of the 
Fermi circle {i.e., nonzero \Ep — Ed\), its period changes as the chemical potential (or the 
doping level) of graphene is changed. 

Although the dangling bond causes the increase in the hopping parameter at the edge, it 
hardly changes the proportion of two backscattering channels, and the electron density profile 
of the armchair edge with the dangling bonds (Fig. 2b) is almost indistinguishable from the 
hydrogen-passivated armchair edge (Fig. 2a). For the pentagonal reconstruction (called a 
5-6 configuration) at the armchair edge (Fig. 2c), both the intervalley and the intravalley 
scatterings take place for the obliquely incident wave because of different boundary matching 
conditions from the ideal edge, which means that the orientation of pseudospin now changes. 
Although this mixing of intervalley and intravalley scattering channels slightly modifies the 
charge modulation pattern, overall features of the electronic density profile in the right panel 
of Fig. 2c looks similar to those of the hydrogen-passivated armchair edge. In particular, we 
note that the node-like structure survives here. The reason is that the small-angle incident 
waves have large weights when forming a standing wave near the edge and are reflected 
dominantly via intervalley scattering for the armchair edge regardless of the details of edge 
termination. 

At the zigzag edge, on the other hand, only the intravalley scattering is allowed which 
shows a standing wave quite similar to the conventional metal surface. The electronic density 
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profile of the zigzag edge is given by 

7T 

= 2 §\ | exp(— ik x x) — exp(ik x x))\ 2 d9 

oc 1 - J (2k R x). (2) 

Here, to reveal the characteristic feature of the scattering unambiguously, the contribution 
of the edge states at the zigzag edge is not included in the electronic density profile. With 
the intravalley scattering at the zigzag edge, the pseudospin changes its orientation. Since 
the scattering is confined to the same valley, the pseudospin follows the change of the crystal 
momentum by scattering. This is the sharply constrasted behaviour of pseudospins at the 
armchair edge discussed above. 

For a low-energy scattering, the diameter of the Fermi circle in graphene is small. Specif- 
ically, the Fermi circle diameter of slightly n-doped epitaxial graphene on the SiC(OOOl) 
surface is about 0.1 A" 1 so that the intravalley scattering makes a long-period charge mod- 
ulation pattern of a characteristic period of ~3 nm. Figure 3a is the real space topographic 
image and Fourier-filtered line profile of the experimental STM data of the armchair edge 
(see the Methods section). By Fourier-filtering, it is possible to suppress other scattering 
effects from the substrate or from some unidentified defects which may exist in graphene. |21j 
The DFT-simulated STM image for the hydrogen-passivated armchair edge in Fig. 3b agrees 
excellently with the experimetal result (Fig. 3a). We have also simulated the armchair edge 
with the 5-6 configuration [T8] and demonstrated that the atomic scale node-like structure in 
the y-axis direction also occurs, as shown in Fig. 3c. These are contrasted to the case of the 
zigzag edge in Fig. 3d where no node-like structure but a periodic modulation is observed. 

Next, we move on to scattering properties of graphene with super-periodic edge structures 
with atomic-scale defects in Fig. 4. In reality, graphene sheets usually have such complex 
edge structures. The super-periodic armchair edges have no other propagating backscatter- 
ing channels than the original intervalley and intravalley ones, if the radius of the Fermi 
surface is sufficiently small. A double-vacancy case and a single-adatom case for different 
supercell periods are considered. The solid (dashed) lines in Fig. 4a and 4b are reflection 
probabilities of the armchair edge with relatively high (low) defect density, and the reflection 
probability is expected to approach asymptotically to the perfect armchair case in the low 
density defect limit. The overall shape of the reflection probability is quite insensitive to 
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the hopping parameters of edge sites and the incident energy. 

When vacancies are introduced in the zigzag edge, on the other hand, the scattering 
channels change drastically depending on the period of the edge. When the period becomes 
three, the intervalley channel opens because both K and K' fall on the T-point in the folded 
Brillouin zone, and it becomes a dominant backscattering channel. If a carbon atom is at- 
tached to the zigzag edge, the scattering becomes rather complicated. Though the reflection 
probability converges to the perfect zigzag edge limit at a low density of adatoms, a high 
density of adatom makes the reflection probability strongly depend on hopping parameters 
and the incident electron energy. It is also found that there exists an abrupt suppression 
of one channel around the incident angle of 90° in the super-periodic zigzag edge, which is 
due to the asymmetric trigonal warping [22] between K and K' valleys. This asymmetry 
renders only the intravalley scattering channel available [23J for a certain k y range for the 
zigzag edge, and that property was exploited for generating valley-polarized currents in the 
former theoretical work. j21] Such a suppression is negligible at the armchair edge because of 
the mirror-symmetric pattern of trigonal warping with respect to the k y -axis. (See the Fig. 
S2 of Supplementary Information for the warping of K and K' valleys for the armchair and 
zigzag edges.) More detailed analysis and understanding of the scattering at realistic edges 
with defects certainly require enormously more work in the future. 

Methods 

Theory We used the density functional theory [25] [26] calculations and the tight-binding 
boundary matching methods to obtain the real-space STM images and reflection probabil- 
ities at graphene edges. The Perdew-Burke-Ernzerhof functional formj[2"T] was adopted in 
the generalized gradient approximation and the ionic potentials were described by projector- 
augmented waves [28] (PAW) implemented in the Vienna Ab initio Simulation Package 
(VASP).[29j The plane wave basis with the kinetic energy cutoff at 300 eV was employed for 
describing wavef unctions, and graphene edges were fully relaxed until the force on each atom 
was within 0.02 eV/A. Tight-binding parameters were obtained from maximally localized 
Wanneir function analysis using wannier90.[3T)] 

Experiment The graphene was grown epitaxially on the silicon face of a highly n-doped 
6H-SiC(0001) by thermal desorption of silicon at high temperature. To obtain high quality 
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graphene and its edge structure, the pressure during the thermal desorption was kept below 
3 x 1CT 9 Torr. The STM images were obtained at 300 K in ultrahigh-vacuum with an 
Omicron instrument. In monolayer graphene, all hexagons were clearly seen and we could 
easily confirm the crystallographic direction. In order to filter out other effects such as the 
substrate structure and the graphene honeycomb lattice structure coexisting in the STM 
image, we applied the Fourier Filtering in the WSxM software. 
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FIG. 1: Scattering channels at the graphene edge. a,b, Schematic Fermi surfaces and pseu- 
dospin fields of the armchair edge (a) and the zigzag edge (b) . The Brillouin zone is transformed to 
a rectangle in each case, for convenience in taking into account the translational symmetry in the 
y-direction alone by the presence of the edge. Green arrows stand for orientations of pseudospin 
fields. c,d, Available scattering channels of a given incident wave for the armchair edge (c) and 
the zigzag edge (d). Red arrows indicate the incident direction and blue arrows the scattering 
direction, respectively. In d, the incident wave is assumed to be .fT-valley polarized. In c, 9 and k 
denote the incident angle and the Bloch wavevector, respectively. The size of the Fermi circle is 
exaggerated for visual clarity. 
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FIG. 2: Reflection probability of the K-valley polarized incident wave and the laterally 
averaged electronic density profile for different edge geometries. a,b,c,d, Reflection prob- 
abilities of the if-valley polarized incident wave (central panels) and laterally averaged electronic 
density profile from both valleys (right panels) at the hydrogen-passivated armchair edge (a), the 
dangling armchair edge (b), the reconstructed armchair edge (c), and the hydrogen-passivated 
zigzag edge (d), respectively. The Bloch cell-periodic part is ignored for clarity. The single hop- 
ping parameter of —2.88 eV is sufficient to describe the hydrogen-terminated armchair graphene 
edge. At the edge of the dangling bond, however, our Wannier function analysis demonstrates 
that hopping parameters are increased by 30% along the edge. At the edge of the pentagonal 
reconstruction structure in c, on the other hand, hopping parameters are reduced by 30 — 50%. 
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FIG. 3: STM topographic image of the graphene edge, a, Experimentally observed STM 
topographic image of the graphene armchair edge and its Fourier-filtered profile along the black 
straight line in the image. b,c,d DFT-simulated STM images for the hydrogen-passivated armchair 
edge (b), the 5-6 configuration edge (c), and the zigzag edge (d). The intervalley scattering 
manifests itself in a clear node-like structure of the charge density in the y-axis direction in b and 
c, and the DFT-simulation in b is in excellent agreement with the experimental image in a near 
the edge. 
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FIG. 4: Reflection probability of the iT-valley polarized incident waves. a,b,c,d, Reflec- 
tion probability of the K- valley polarized incident wave of double vacancies at the armchair edge 
(a), adatoms at the armchair edge (b), single vacancies at the zigzag edge (c), and adatoms at the 
zigzag edge (d). Fractional numbers indicate the number of defects per unitcell along the edge. 
For the super-periodic zigzag edge with defects, the intervalley scattering channel is open when 
the period is a multiple of three, because K and K' points coincide with the T-point when the 
zigzag Brillouin zone is folded three times. Therefore, both intervalley and intravalley scatterings 
are allowed for the 3m (to: an integer) period. Since only the intravalley scattering is allowed at 
the (3m+l) and (3to+2) period edges, the plots of reflection probabilities for these cases are not 
presented. 
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